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Abstract 

Working with a magnetic field periodic along Oz and decaying in 
time, we deal with the Dirac-type equation characterizing the fermions 
evolving in magnetar's crust. For ultra-relativistic particles, one can 
employ the perturbative approach, to compute the conserved current 
density components. If the magnetic field is frozen and the magnetar is 
treated as a stationary object, the fermion's wave function is expressed 
in terms of the Heun's Confluent functions. Finally, we are extending 
some previous investigations on the linearly independent fermionic 
modes solutions to the Mathieu's equation and we discuss the energy 
spectrum and the Mathieu Characteristic Exponent. 
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1 Introduction 



Many years after Duncan and Thompson, [Duncan and Thompson, 1992], 
introduced the notion of magnetars, for almost non-rotating neutron stars, 
whose magnetic field strength was estimated to be about 10^ to 10^ larger 
than the one of a neutron star, this class of astrophysical objects has become 
an active of research. 

With magnetic fields greater than the critical induction at which the 
cyclotron energy of an electron equals its rest mass energy, they are the only 
source of the brightest cosmic events originating outside the solar system, 
known as giant flares. 

The magnetars exotic properties have led to many open questions, espe- 
cially related to the configuration of the magnetic field inside and to their 
internal structure [Sinha and Mukhopadhyay, 2010; Orsaria et al., 2011]. 

Today, it is no doubt that the crust/core couphng [Glampedakis and 
Andersson, 2006], the elastic properties of the outer part of a neutron star 
[Pethick and Potekhin, 1998], the magnetic field geometry and the effect 
of the strong magnetic pressure on the tension in the star crust, which are 
all alterating the detected modes frequencies, are strongly depending on the 
crustal composition [Piro and Bildsten, 2006] , and hence on the equation of 
state (EOS) [Lattimer and Prakash, 2004]. 

Theoretically, by solving the Einstein equations for a general time-invariant 
spherically symmetric metric, one comes to the celebrated Tolman-Oppen- 
heimer-Volkoff (TOV) equation which is expressing the dependence of the 
star matter density and pressure on its mass and radius. Reversely, maximum 
masses and minimum periods of the neutron stars allow us to set constraints 
on the EOS. 

The present paper is following some previous investigations on spinless 
particles, described by the Klein-Gordon equation, moving in strong mag- 
netic induction periodic along Oz [Dariescu et al, 2011a; Dariescu and Dari- 
escu, 2011b]. 

Since magnetars are being composed of different types of fields which all 
have deep implications on their parameters [Crawford and Kazanas, 2009], 
the present work is a general analysis of fermions evolving in the mag- 
netic field configuration proposed by Wareing and HoUerbach, in their three- 
dimensional pseudo-spectral numerical MHD code [Wareing and HoUerbach, 
2006; Wareing and HoUerbach, 2009]. 
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2 Ultra-relativistic Particles in the Crust 



Generally, the covariant theory starts from an effective Lagrangian of the 
nucleons and as many meson fields as needed to produce the basic nuclear 
matter. 

By neglecting the self- interaction terms, one comes to the decoupled Dirac 
and Klein-Gordon-type equations (in natural units h = 1 = c) 



and 



YDi + mo * = , 
$ = 0, 



where the covariant derivatives, 

Di = di- iqAi , 



(1) 
(2) 



include the interaction with the external fields, Ey(z,t) and Br^(z,t), gener- 
ated by the Ay component of the four-potential. 

By expressing the bi-spinor "if in the equation (1) as 



mo 



X: 



(3) 



and using the Dirac representation for the {Y}ti inatrices 
with 



a'- 



af" 
a^" 



where X is the unit 2x2 matrix and o"^ are the usual Pauli matrices, the 
equation (1) leads to [Dariescu and Dariescu, 2012] [M. A. Dariescu and C. 
Dariescu, submitted to Astroparticle Physics] 

This contains the matrix Fi = ^7^7^, whose eigenvectors 

1 

±1 
1 

±1 



(4) 



^1,2 = - 



(5) 
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arc corresponding to the eigenvalues Ai^2 = =tl- 

The explicit form of the equation (4) is depending on the configuration of 
the background fields in the magnetar's crust. This is a very active field of 
research and the expression of the strong magnetic induction is not unique. 
One may check easily that the following expressions 



a 



bo sin (kz) exp 



Bx — bo cos (kz) exp 



t 

a 

— t 
a 



(6) 



sustained by only one component of the four-potential, namely 

bn «^ 

Ay^- — sin {kz) e""* with a = — (7) 

K/ (7 

are solutions of the Maxwell's Equations 

dB ^ - ^ ^ dE 

VxE^-—, V-E^O, V-B^O, VxB^j + — . 

at at 

For a magnetic field frozen in the crust, the time-independent part of the 
magnetic induction in (6) is of the form proposed by Wareing and Hollerbach 
[Wareing and Hollerbach, 2006; Wareing and Hollerbach, 2009], i.e. B^ = 
BQh{z), where h{z) is a periodic function. 
With the standard variables separation 



Xi,2 = ViA^, t) exp [ip± ■ x±] wi,2 , 



(8) 



where p± ■ x± = p^x + Pytj and the index "1" or "2" stands for Ai = +1 and 
A2 = —1, the equation (4) leads to 



dz^ 



dt' 



p'i + ml + 2py— sm{Kz)e~"^ + 



K 



'qbo^ 



K 



Biv?[Kz)e 



-2at 



±qbo cos{Kz)e " 7] 



0. (9) 



Analytically solving of the above equation in the general case when all 
the terms are taken into account is a difficult task and therefore we need 
some simplifying assumptions. 
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For ultra-relativistic particles with the momentum along the electric field, 
Py ^ qbo/K, so that {py — qAy)^ ft! p^, the equation (9) simplifies to 



2 I 2 



and one may use the Laplace-Fourier expansion 



T) ± qbocos{Kz)e " rj — 



-{na+iu!)t 



(10) 



(11) 



n=0 



to get the following recurrent differential system for the amplitude functions 



+ 



P 



— 2niauj — r?o? 



Ipn ± qbo COs{kz) = , 



(12) 



where = cu^ — — ml. 

For n — 1, the relation (12) becomes the inhomogeneous second order 
linear differential equation with constant coefficients, 



1^ 



+ 



p^ — 2iauj — ip ± qbo cos{kz) tpo — , 



where 



is satisfying the initial differential equation 



dz^ 



Using in (13) the standard decomposition 

we are identifying the coefficients, 

r qbo >^ + 2/«p + — liaijj ^ qbo np — iau 

A — ±J\f --— -— r — -r ~ ±7V 



(13) 
(14) 

(15) 
(16) 



2 (k"^ + 2kp + a'^y + A(aujY 

2 («2-2«p + Q;2)2 + 4(Q;a;)2~ ' 
5 



4 {Kpf + {aujf' 
qbo —Kp — iauj 



4 {kpY + (auj) 



2' 



and obtain the following mode expressions 

V'l,2(^) ^ 



qbo Kpsm{Kz) — acu cos{kz) 
1 (fi;p)2 + (aa;)2 



2 ^J{kpY + (aa;)^ 



sin - ^) e*^^ , 



;i7) 



where 



tan 6' 



up 



is a phase factor in the wave functions of chiral ultra-relativistic fermions of 
energy cu, evolving in the magnetic field = 6oe "*cos (kz). 

Going all the way back, the two linearly independent modes of positive 
energy, solutions of (1), are 



Px — iPy + p — uj + ia — mo — ih{z) 
Px + iPy — p — oj + iot — mo + ih{z) 
—Px + ipy — p + UJ — ia — mo + ih{z) 
—Px — iPy + p + UJ — ia — mo — ih{z) 



i/jie 



[i{p±-x±-ujt)] -at 



(18) 



and 



*2 = 



—Px + ipy + p — UJ + ia — mo — ih{z) 

Px + iPy + p + UJ — ia — mo — ih{z) 

Px — iPy — p + UJ — ia — tuq + ih{z) 

—Px — iPy — p — UJ + ia — mo + ih{z) 

where we have introduced the notation 

kpcos{kz) + Q;a;sin(«;2;) 



^2 6 



[i{p±-x_L-u>t)] ^-at 



(19) 



h{z) = K 



Kpsm.{Kz) — auicos{Kz) 



K cot {kz — 9) . 
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For the superposition 



Px - iPy 

-p — uj + ia + ih{z) 

-Px + iPy 



^^^[iip±-S±-u^t)] ^-at ^ ^20) 



p + CO — ia — ih{z) 

the current density components defined by 

j = , / = g = g*t^, 

have the exphcit expressions 

= q[{uj+p)p^ + {a + h)py\\ipi\'^e-'^'^^ 
jy = q[{uj+p)py-{a + h)p^]\tlJi\^e-^''' 



(21) 



= ^[2p'' + 2ujp + ml + {a + hf 
Q = ^\2uj^ + 2ujp-ml + {a + hf 



(22) 
resonance 
(23) 



which are fulfiUing the conservation law dif = 0, for a = and the 
condition 

^ + 2up + mo . 

For q; 7^ 0, at 2; = L^, once the longitudinal momentum characterizin 
the particles confined in the crust is equal to 

p = — cot [kLz) , 

K 

i.e. kLz — 9, the wave function absolute value and consequently jx and j 
are all vanishing, while the other current's components, at 2; = L^, i.e. 

Uz = L.) = Qiz -Lz)-l \^ff ' ^ sin^ i^L,) e"^"* , 

are satisfying dzjz + dtQ = 0. 
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3 HeunC Solution within Ohmic Timescale 



In the opposite case, corresponding to Py <C qbo/n and the time variable 
less than the characteristic time, so that e~"* ~ 1, which is a more realistic 
situation for a magnetar characterized by a strong magnetic induction, bo ~ 
10^^ G, the four-potential component (7) can be taken as 



Ay = sm KZ . 

K 

With the function r]{z,t), in the general equation (9), written as 

r)(z,t) ^v{z)e-^\ 



(24) 



the chiral components </?i,2(2;), corresponding to A = ±1, are solutions of the 
equations 

where ( = kz and qbo = b. 

In terms of the new variable s = cos(, that maps the ( space to s e 
(-1, 1), with 

b 



<^i,2(s) = exp 



in (25), we come to the following differential equation for u, 



+ 



ds 



p 







(26) 



(27) 



Up to normalization constants, the solutions of (27) are the HeunC confluent 
functions [Arscott, 1995] 

r 7i 1 1 7i 5^ r>2 

iJ'i,2{s) e < HeunC 



1 1 „63 p' b s 
±4—, --, --, t2^, -:-^± 



2' 2' ' 8 ^2 ^2' 2 



2(1 + s) HeunC 



2' ' 8 «;2 
'b I 1 63 p2 h s + l 
2' ~2' ^2^' 8 ~ ^ ^ ^' 



(28) 

with the boundary conditions u{s = ±1) = and u'{s — ±1) = finite, so 
that the wave functions (26) have the explicit expressions 



exp 



b 

±— cosC 



Ml,2 [cOS^ ^ 



(29) 



where Ui^2 a-re given in (28). 

Now, one may follow the procedure developed in [Gurappa and Panigrahi, 
2004], for a differential equation which can be cast in the form 

[F(D) + P(x, d/dx)] u{x) = , 

where D = x^, F(D) — Y.n'^nD'^ is a diagonal operator in the space of 
monomials and P{x, d/dx) is an arbitrary polynomial function of x and d/dx. 
The necessary condition for a polynomial form of u is 



leading, for 



F(D)x" = , 



in (27), to the following longitudinal momentum quantization law 



p 



2 2^b' 

n K. + — , 



or, for 6 ^ K, to the familiar resonance relation p = hk. 

Even though we have derived an analytic solution to the Dirac equation 
(1) for the fermions minimally coupled to the electromagnetic field generated 
by the potential component (24), we must signal the fact that there are un- 
solved problems regarding the normalization of the Heun functions and their 
behavior for some domains of the complex variable. In terms of available soft, 
their treatment can be done only with MAPLE routines which are working 
well in particular cases, but are very sensitive to the parameters, breaking 
down for values which might be of interest. 

Finally, if we keep only the first order contribution in b/n^, the equation 
(25) becomes the well known Mathieu's equation [Gradshteyn and Ryzhik, 
1965] 

d'^v 

—^ + [a-2/3 cos(2«;)] y 



dw"^ 



0, 



(30) 



with w — Kz/2 and 



2b 



whose solutions are the so-called Mathieu's even and odd functions C and 5, 



C 



4p^ 2b KZ 

- 2 ' ~^~Z2 ' ~0 



K 



s 



4p^ 2b KZ 
- 2 ' "F^i ' "o" 



K 



(31) 
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These can be written as 

where / is a periodic function and the so-called Mathieu Characteristic Ex- 
ponent (MCE), 7, is strongly related to the particles momentum along Oz 
and to the external magnetic field configuration. 

For a » ^, the solution of (31) is stable and the energy spectrum remains 
continuous almost everywhere. In our case, this means ultra relativistic par- 
ticles with high momentum along Oz, i.e. 

22 2 22 24 -^-^ HqbQC 

P C ^ CO - p^C — TJIqC > — - — , 



where hqboc^ 4 {MeVy, for bo ~ 10^^ G, which is typical for the magnetic 
induction in magnetar's crust. 

Reversely, for fcrmions characterized by ^ ^ a, the periodic term in (31) 
can no longer be treated as a perturbation, and one gets discrete energy 
eigenvalues, separated by large instability regions, which become wider as /3 
increases. The asymptotic expansion [Beckert et al, 1979] 

a^-2l3 + An^, (32) 
leads to the following energy quantized levels 



u!^ Ki p']_ + ml + nK^2qbQ ± qbo . (33) 

Moreover, for particles whose longitudinal momentum is close to k/2, in 
the range 

Kb K b 

<p< --\ , (34) 

2 2/t 2 2k ' ^ ' 

the imaginary part of MCE comes into play, leading to exponentially growing 

wave functions. 

This can be seen by writing down the solution of the general Mathieu's 
equation (31) as 

n 

where the C— coefficients are satisfying the relation 

aC2n - (7 + 2n)^C2n " /3C'2(n-l) " /3C'2(n+l) = , 
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or, m our case, 

C2n + i2n (^2(71+1) + C'2(n-1) = , with ^2n = ± 



2h 



4p2 - (-y + 2n)2«;2 



By imposing the vanishing of the determinant associated to the above ho- 
mogeneous system of unknown C_2(n+i), • • • , C'2(n+i), we get the following 
relation for the MCE 



cos(7r7) = 1 - 2A(0) sin^ 



P 

TT — 



(35) 



where A(0) is the value determinant if we set 7 = and p 7^ nn. Generally, 
7 is real or imaginary, depending on the model parameters. 
For p in the range (35), i.e. 



TT 

2 



P 71 
< TT- < - 
K 2 



1 + ^ 



as soon as 7 becomes imaginary and |cosh(77r)| > 1, the wave function is not 
bounded on the real axis [Grib et ai, 1994], while for the resonance condition 
p = nn, the relation (36) is replaced with [Coisson et ai, 2009] 

cos(7r7) = 2A(1) - 1 . 

For bosons, we have come to similar results [Dariescu and Dariescu, 2011], 
they being described by the wave function 



i( P±-x±—Qt 



where 



c 



KZ 



2k^ ' 4«;4 ' 



KZ 



and by the energy quantization law 

Q2 + + 2n qbo . 



(36) 
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4 Conclusions 



As in our previous works, [Dariescu et ai, 2011a; Dariescu and Dariescu, 
2011b], we have considered that the periodic magnetic field proposed by 
Wareing and HoUerbach [Wareing and HoUerbach, 2006] is hkely to exist in 
magnetar's crust and we have analyzed how the time and spatial distribTition 
of the external fields (6) have an infiuence on the wave function, solutions to 
the corresponding Dirac equation. 

In a perturbative approach, for ultra-relativistic particles, we have derived 
the current density components and the charge density (22). These can 
be used for studying the back-reaction of the new fields generated by the 
Maxwell equations. 

Within the Ohmic timescale, when the magnetic field is frozen in the 
crust and the magnetar is treated as a stationary object, the fermions wave 
functions can be written in terms of the Heun's Confiuent functions. 

Today, different types of Heun's functions are considered as successors 
of the hypergeometric functions, with a wide application in modern physics 
and new algorithms have been created, for finding solutions to systems of 
nonhnear transcendental equations [Fiziev and Staicova, 2012]. For example, 
the confluent Heun functions have been worked out in the context of the 
quasinormal modes for nonrotating black holes [Fiziev and Staicova, 2011] or 
as solutions to Schrodinger equation, for different rational potential functions, 
in thick braneworlds [Cunha and Christiansen, 2011]. 

The flnal part of the present paper, which extends some previous works 
[Dariescu et al, 2011a; Dariescu and Dariescu, 2011b], is leading to conclu- 
sions in agreement with the linear stability approach, belonging to Rhein- 
hardt and Geppert [Rheinhardt and Geppert, 2001], who have proved that 
the transfer of magnetic energy from the background field to small-scale 
modes may produce Hall instabilities, at some ranges of wavenumber values. 

As we have mentioned previously, the configuration of the background 
magnetic field in the magnetar's crust has been extensively investigated and 
is not unique. For example, for a magnetic field oriented along Ox and 
everywhere parallel to the surface of a slab with a finite 2;— thickness, the 
induction has been written as = B^f^z), where f{z) is vanishing for 
z > L. Guided by the boundary conditions and working with dimensionless 
variable z/L, the authors of [Rheinhardt and Geppert, 2002] have proposed 
the form / = (l — j^). With this choice, the ultra-relativistic particles 
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would be described by the differential equation 



1^ 



1 - 



L 



^ = 0, (37) 



where h = qBo and = oj"^ — p']_ — m^. This is a particular case of the 
analysis developed in the present work since one can easily come to the 
above equation by performing, in (10) written for a ?s and 77 = ■0(2;)e~"^*, 
the series expansion 



cos(«;2;) fa 1 — 



2 ' 

valid for 2; <C L and by making the identification \/2/k L. 

As a general remark, we would like to mention that the present work 
could be extended in several directions. 

Thus, because besides leptons, in dense nuclear matter, the nucleons inter- 
act among each other through three meson fields: the isoscalar-scalar meson 
(T, isoscalar vector meson u and isovector-vector meson p, one may assume 
that currents and charge density produced by the single-particle spinors, as 
the ones in (22), are actually acting as sources for the Klein- Gordon Gordon 
equations for the time- and space-like meson fields 



ml) 00 = 94>Q' (38) 



where 4> — {uj , p) and 0o — {^0 ; Po}) g<t> ^'^'^ the respective coupling con- 
stants, are the meson masses and (see (21)) 



1=1 i=l 



for a number of protons Z. 

In the Hartree approximation and considering only static configurations, 
the system (39) is producing new meson fields and one should look for a 
closed set of solutions which usually is found itcratively, until convergence 
results. For a detailed study of the effect of strong magnetic field on nuclei 
in the crust of magnetars and its consequences on the physical measurable 
quantities, we recommend [Artega et ai, 2011]. 

Last but not least, the rclativistic particles energy spectrum can be em- 
ployed in dealing with the Equation of State and its implications on the 



13 



detected modes frequencies. This possible extension is very challenging due 
to the intense activity on several types of EOS, within different theoreti- 
cal techniques and for different compositions of the neutron star, which can 
account for the observed frequencies and magnetic field strengths. 
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